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Abstract 

We refine and update the metastability constraint on the Standard Model top and 
Higgs masses, by analytically including gravitational corrections to the vacuum decay 
rate. Present best-fit ranges of the top and Higgs masses mostly lie in the narrow 
metastable region. Furthermore, we show that the SM potential can be fine-tuned 
in order to be made suitable for inflation. However, SM inflation results in a power 
spectrum of cosmological perturbations not consistent with observations. 



1 Introduction 

Assuming that the Standard Model (SM) holds up to some high energy scale close to Mp\ = 1.22 10 19 x 
GeV, present data suggest a light Higgs mass, ~ (115 — 150) GeV. If the Higgs is so light, radiative 
corrections induced by the top Yukawa coupling can destabilize the Higgs potential and the electroweak 
vacuum becomes a false vacuum, which sooner or later decays [1, 2, 3, 4, 5, 6]. Demanding that the 
SM vacuum be sufficiently long lived with respect to the age of the universe implies a bound on the 
Higgs and top masses [2, 3, 4, 5]. 

In section 2 we recall the peculiarities of vacuum decay within the SM relevant for our later 
inclusion of gravity, which was neglected in previous analyses. In section 3 we show how gravitational 
corrections to the vacuum decay rate [7] can be computed by making a perturbative expansion in the 
Newton constant, and we obtain the analytic result for SM vacuum decay. 

In section 4 we show that for fine-tuned values of the Higgs and top masses (that lie within 
the experimentally allowed range), the SM potential can be suitable for inflation. However, the 
corresponding power spectrum of anisotropies is larger than the observed one. 



2 Vacuum decay within the Standard Model 

We recall vacuum decay within the Standard Model without gravity, and its peculiarities relevant for 
our later inclusion of gravity. The SM contains one complex scalar doublet H, 



H 



(h + iri)/\ 



(1) 



with tree-level potential 



V = m 2 \H\ 2 + A|F| 4 = -m 2 h 2 + \\h 4 + ... (2) 
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where the dots stand for terms that vanish when the Goldstones r], x are set to zero. With this 
normalization, v = (G F V2)~ 1/2 = 246.2 GeV, and the mass of the single physical degree of freedom 
h is rn? h = V"(K)\ h=v = 2\v 2 . As is well known, for h S> v the quantum corrections to V(h) can 
be reabsorbed in the running coupling \(fl), renormalized at a scale /Z ~ h. To a good accuracy, 
V(h 3> v) ~ A(/i)/i 4 /4 and the instability occurs if A becomes negative for some value of h. For the 
values of rrih compatible with data this occurs at scales larger than 10 5 GeV, suggesting that we can 
compute vacuum decay neglecting the quadratic term m 2 h 2 /2. 

The bounce [1] is a solution h(r) of the Euclidean equations of motion that depends only on the 
radial coordinate r 2 = x^x^: 

-d^h + VXh) = -^-lf + V>(h)=Q, (3) 

and satisfies the boundary conditions 

ti(0) = , h(oo) = v -> . (4) 

We can perform a tree-level computation of the tunnelling rate with a negative A < renormalized at 
some arbitrary scale fx. In this approximation, the tree-level bounce ho{r) can be found analytically 
and depends on an arbitrary scale R: 

Mr) = ■ 5 oN = ~ • (5) 

At first sight, computing the decay rate among two vacua in the approximation V(h) = A/i 4 /4 where 
no vacuum exists may appear rather odd. However [4], the presence of a potential barrier around 
the false vacuum h ~ is not necessary, since in quantum field theory the bounce is not a constant 
field configuration, and the energy in its gradient effectively acts as a potential barrier. Furthermore, 
the decay rate does not depend on the unknown physics that eventually stabilizes the true vacuum 
at h ~ Mpi, if the bounce has size R S> 1/Mp\: once a tunneling bubble appears, the instability due 
to V'(h(0)) 7^ brings h down to the true minimum with unit probability. Formally, by performing 
the analytic continuation from Euclidean r 2 = x 2 + t 2 to Minkowskian r 2 = x 2 — t 2 space-time, the 
evolution inside the bubble is described by eq. (5) with r 2 < 0: ho reaches a singularity at r 2 = — R 2 . 
Indeed our potential is unbounded from below. 1 In general, what happens inside the bubble does not 
affect the tunneling rate nor the growth of the bubble: being an 0(4)-invariant configuration (i.e. the 
bounce depends only on r), its walls expand at the speed of light, so that what happens in the interior 
cannot causally affect the exterior. 

The arbitrary parameter R appears in the expression of the SM bounce ho(r) because in our 
approximations the tree level SM potential is scale-invariant: at this level, there is an infinite set of 
bounces of varying size R, all with the same action So [ho] . 

Quantum corrections are the dominant effect that breaks scale invariance, and have been computed 
in [4]. At one-loop order, the tunnelling probability in the universe space-time volume Vjj is then given 
by 

p = m^[p(R)], p(R) = ^ e - s , (6) 



In the usual case, with a potential with two minima, the bounce can be computed only numerically. The analytic 
continuation can be done by switching r — * ir in eq. (3) at r < 0, and solving numerically. The qualitative behavior of 
the solution can be understood by noticing that this operation is equivalent to flipping the sign of V: h oscillates around 
the true minimum, reaching it at r — » — 00, i.e. for asymptotically large times inside the expanding bubble. 



2 



where S = Sq + ASi_i oop is the one-loop action: since the bounce is not a static configuration, 
corrections to both the potential part, as well as to the kinetic part of the action, must be taken into 
account [4]. To find the bounce configuration that extremizes S, it is enough to evaluate it along 
the family of tree-level bounces, ho in eq. (5), and minimize with respect to R. The result is roughly 
S ~ 8-7r 2 /3|A(/2 = 1/R)\, i.e. one-loop corrections remove the tree-level ambiguity on the RGE scale p 
by fixing it to be the scale 1/R of the bounce. Since within the SM A(/2) happens to run reaching a 
minimal value at p ~ 10 16_17 GeV, tunneling is dominated by bounces with this size. A posteriori, 
this justifies having neglected the SM mass term, that gives a correction AS ~ (mR) 2 <C 1 to the 
bounce action, and suggests that gravity should be taken into account perturbatively. 



3 Vacuum decay with gravity 

We now extend the previous computation taking into account gravity [7]. In our case this is a 
potentially relevant effect, since gravity breaks scale-invariance and 1/R is just somewhat smaller 
than the Planck scale. One might worry that gravity can have dramatic effects, and that the decay 
rate starts to depend on the unknown depth V m i n of the true minimum of the SM potential. 2 This is 
not the case. Since the exterior geometry is the flat Minkowski space, the generic argument given in 
the non-gravitational case still holds: the bubble is an 0(4)-invariant solution and its walls expands 
at the speed of light, irrespectively of what happens inside. 3 

We recall from [7] the basic formalism needed for a quantitative analysis. We assume an Euclidean 
spherically symmetric geometry, ds 2 = dr 2 + p(r) 2 d£l 2 , where df2 is the volume element of the unit 
3-sphere. The Einstein-Higgs action 



S 



d 4 x^/g 



(dM&*h) 



+ V(h) 



2 k 



(8) 



where re = 8vrG and G = l/M^ with M P i = 1.22 10 19 GeV, simplifies to 



S = 2ir 2 



dr 



b' 2 3 
P (-7T + V) + -0 p +pp 

2 K 



p) 



(9) 



where ' denotes d/dr. The equations of motion are 



h" + ?> P -ti 
P 



dV 
~d~h : 



'2 



P 



re 2f h' 2 
1 + -p 2 ( — 
3 P v 2 



V). 



(10) 



We can analytically include the effect of gravity, assuming RMp\ S> 1, by performing a leading- 
order expansion in the gravitational coupling re: 



h(r) = h (r) + re/ii(r) + C(re 2 ), p(r) = r + repi(r) + C(re 2 ). 



(11) 



The action is 

S = Sq + 67r 2 re 



dr 



r 2 pi 



h' 2 



+ V(h ) + (rp' 2 + 2p lP \ + 2 Pl rp'[) 



+ 0(re 2 ). 



This is what one would naively guess from the results of [7] for the bounce action: 

Swith gravity ^ ^without gravity 



(12) 



(7) 



i.e. the bubble does not exist if the true minimum has a large negative cosmological constant, e.g. Vmm ~ — Mpi. 
However, eq. (7) holds within the thin- wall approximation [7], not applicable when V m i n is large and negative, and not 
applicable to the SM case we are interested in. 

3 It is only an observer inside the bubble that experiences a large negative cosmological constant and consequently a 
contraction down to a big-crunch singularity [7], instead than an expanding bubble. 
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Figure 1: Probability p{R) that the SM vacuum decayed so far for rrih = 115 GeV, nit = 174.4 GeV, 
otz(Mz) = 0.118, due to bounces with size R, without including gravitational effects (dashed 
curve [4]) and including gravitational effects (continuous line). The correction is relevant only at 
1/R > 10 17 GeV. Uncertainties due to higher-order corrections are not shown. 



We have taken into account that many terms in the expansion vanish either because the integrand is 
a total derivative (e.g. the negative power 1/ k in eq. (8) is just apparent) or thanks to the equations 
of motion. Indeed h\ does not appear in eq. (12) because we are functionally expanding around the 
extremum ho of the non-gravitational action, so that the first functional derivatives vanish thanks to 
the equations of motion. So, we only need to compute p±: its equation of motion is 

A = \r 2 (f " V(Ao)) • (13) 

Inserting it into eq. (12) completes the computation of gravitational corrections to leading-order in 
k. We notice that the first term in eq. (12), which is linear in p±, contributes —2 times the last 
purely gravitational term in eq. (12), which is quadratic in p\. This happens because S must have 
an extremum at c = 1 under the variation p\{r) — > cp\{r). The discussion is so far general, and by 
choosing toy potentials we verified that eq. (12) agrees with the full numerical result. 



4 Vacuum decay with gravity in the Standard Model 

Going to the SM case, using the analytic expression of eq. (5) for the bounce ho, we can perform all 
integrations analytically finding 

87T 2 _ 2567T 3 

3\X\ + A6l - lo ° P + A, Wit y , ^gravity - 45 ( RMp ^2 ^4) 

where AS^^ty is the gravitational correction and A5i_i oop the one-loop correction, given in eq. (3.3) 
of [4]. Eq. (6) gives the tunneling probability p(R). 

Fig. 1 shows an example of the relevance of gravitational corrections. We checked that the leading- 
order approximation agrees with the result of a full numerical computation: eq. (14) correctly approx- 
imates the action of the true bounce, and the true bounce h{r) is correctly approximated by ho{r) 
with the value of R that minimizes S. 4 



Here we comment about the comparison between the analytic result in eq. (14) and the full numerical computation. 
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Higgs mass m h in GeV 



Figure 2: Metastability region of the Standard Model vacuum in the (mh,mt) plane, for a s (mz) = 
0.118 (solid curves). Dashed and dot-dashed curves are obtained for a s (mz) = 0.118 ± 0.002. The 
shaded half-ellipses indicates the experimental range for rat and at 68% and 90% confidence level. 
Sub-leading effects could shift the bounds by ±2 GeV in mt- 



Fig. 2 shows the regions in the (m/,,mi) plane where the SM vacuum is stable, meta-stable or too 
unstable. Gravitational corrections only induce a minor shift on the 'instability' border, less relevant 
than present experimental and theoretical (higher-order) uncertainties. The ellipses truncated at 
mh = 115 GeV are the best-fit values for the top and Higgs masses, from our up-to-date global fit of 
precision data, that includes the latest direct measurement of the top mass, mt = (170. 9± 1.8) GeV [8]. 
Present data and computations indicate that we do not live in the unstable region (such that the SM 
can be valid up to the Planck scale), but increased accuracy is needed to determine if we live in the 
stable or in the small meta-stable region. 

Adding to the SM action possible dimension-6 non-renormalizable operators suppressed by the 
Planck scale would give similar corrections to the bounce action. In particular, adding to the SM 
Lagrangian the operators 

AJ% = {-mlimH\ 2 + + c 2 \H\ 2 \D tl H\^j , (15) 

where £ and c\^ are unknown dimensionless coefficients, gives the following correction 

AS '^» - 15(//p,W ( 128 ^ + fXf +4C2 )' <16) 

which can be comparable to the model-independent gravitational effect computed in eq. (14). 

With a typical potential this is a straightforward procedure: the bounce is determined numerically as a compromise 
between classical solutions which under-shoot and over-shoot the true bounce at large r. With a potential close to ft 4 , 
finding the bounce numerically is more involved: with this potential classical solutions necessarily go to zero at large r; 
however, they generically oscillate to zero as 1/r giving a divergent action. The special feature of fto(r) is that it vanishes 
as 1/r 2 giving a finite action. The true bounce should maintain this behavior. In practice, this is achieved imposing a 
vanishing difference between fto(r) and the numerical bounce. The advantage of our analytic approximation based on 
the set of candidate bounces fto(r) with different values of R is that ill-behaved never enter the computation. 
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Figure 3: Examples of fine-tuned SM potentials that might allow inflation. The right handed axis 
shows the value of the slow-roll parameter e that would give the observed amount of anisotropies. 

The values of the coefficients £ and change under field redefinitions and only their linear 
combination entering (16) is physical. Indeed, under H — ► H(l + a\H\ 2 /Mp-^) we have 5 8c\ = 24Aa, 
5c2 = 6a and 5£ = 0; this transformation can be used to set C2 — > 0. Under the Weyl transformation 
of the metric g^ v — > 0^(1 + a|ff| 2 /Mpj) we have 5£ = a/Wtr, 5c\ = 12aX, 5c2 = a; this transformation 
can be used to set £ — > 0. Both these field redefinitions leave A5g ravity invariant. 

To estimate the magnitude of ASg ravity we can thus restrict the attention to only one of the 
three operators in Ajzfg (we choose the \H\ e term), and estimate its coupling using naive dimensional 
analysis. At one loop, graviton exchanges generate the \H\ e operator with c\ ~ g^y as well as the 
A|-ff| 4 operator with coefficient A ~ g^/ir 2 . Here g s is an unknown coefficient which determines if 
quantum gravity is weakly or strongly coupled, with strong coupling corresponding to g s ~ 7r 2 . One 
might therefore argue that c\ ~ Xir, which implies ASg ravit 

rNw ' ^ ^gravity • 

5 Inflation within the Standard Model? 

For nit ~ 173 GeV and nih ~ 130 GeV (i.e. within the experimentally allowed region) both the quartic 
Higgs coupling A and its /3-function happen to vanish, at some RGE scale around Mp\. Is this just a 
coincidence, or this boundary condition carries some message? Some speculations about this fact have 
been presented in [9]. Here we explore a different aspect, namely a possible connection with inflation. 

The quasi- vanishing of both A and /5(A) allows to have a quasi-flat Higgs potential at h ~ Mpi, 
suitable for inflation. Indeed, we can approximate the RGE running of A as 

A(/x ~ h ) ~ A min + In 2 ^ (17) 

around the special value ho where A reaches its minimal value A mui . The constant 7 is related to 
/3(/3(A)) and has the numerical value 7 ~ 0.6 within the SM. The first and second derivatives of the 

5 We do not distinguish between | D^H] 2 and |iif| 2 |7J M _ff| 2 since these operators coincide on the configurations 
we are interested in. 
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SM potential V ~ \{h)h A /4 vanish at h = h* = hoe x / 4 if A m i n = 7/40967T 4 , such that the slow-roll 
parameters 



M|j fV'y _M* i V" 



£ 16vr V V J ' 11 8vr V 

vanish, allowing for inflation. 

The lack of convincing natural models for inflation might indicate that it happens when scalar 
fields, fluctuating along some vast 'landscape' potential generically unsuitable for inflation, encounter 
a small portion of the potential which accidentally is flat enough. This is what might happen within 
the SM. This potential is illustrated in fig. 3, where we do not show the uncertainty due to higher- 
order corrections, which effectively amounts to a ±2 GeV uncertainty in mj. Can this SM potential 
be responsible for inflation and the generation of anisotropies Sp/p? The answer is: not both. The 
basic problem is that the requirement of having enough e-folds of inflation, 

N = dh/M^^ (i8) 

can be met with a small enough e, but this conflicts with the requirement that quantum fluctuation 
of the Higgs inflaton should also generate the observed power spectrum of anisotropies, 5p/p ~ 10~ 5 , 
i.e. 

— » (0.0054M P1 ) 4 . (19) 

Indeed the height V of the SM potential in its flat region is predicted and cannot be arbitrarily adjusted 
to be as low as needed. This result can be understood by either doing explicit computations with the 
approximated potential \{h)h /4, or by looking at the sample SM potentials plotted in fig. 3. For a 
top mass within the observed range, the plateau is at values of h and V 1 ^ which are are somewhat 
below the Planck scale, but Sp/p at iV ~ 60 comes out larger than the observed value. Successful 
inflation and successful generation of anisotropies would be obtained if for some unknown reason the 
potential would remain flat from h ~ h* up to h ~ Mp\. 



6 Conclusions 

In this paper we have refined and updated the metastability constraint on the Higgs mass, assuming 
the validity of the Standard Model up to the highest possible energy scale, A ~ Mp\. In particular, we 
have taken into account gravitational corrections, which were neglected in previous analyses. These 
corrections turn out to be small and calculable in the phenomenologically interesting region of rrih 
close to its experimental lower bound. The updated constraints in the (rrih,mt) plane are reported in 
fig. 2. Among all possible values, the Higgs mass seems to lie in the narrow region which allows the 
SM to be a consistent theory up to very high energy scales, with a perturbative coupling and a stable 
or sufficient long-lived vacuum. Fig. 4 illustrates the constraints on the Higgs mass as function of A, 
and shows that the (meta)stability constraints do not depend on A when it is around the Planck scale. 

We have also shown that the SM potential can be fine-tuned in order to be made suitable for 
inflation. However, the resulting power spectrum of anisotropies is larger than the observed one. 

Acknowledgements We thank Paolo Creminelli and Enrico Trincherini for useful discussions. 
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